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Abstract 

The following comparison rules for the discrete spectrum of the position-dependent mass 
(PDM) Schrodinger equation are established, (i) If a constant mass mo and a PDM m(x) are 
ordered everywhere, that is either mo < m(x) or mo > m(x), then the corresponding eigenvalues 
of the constant-mass Hamiltonian and of the PDM Hamiltonian with the same potential and 
the BenDaniel-Duke ambiguity parameters are ordered, (ii) The corresponding eigenvalues of 
PDM Hamiltonians with the different sets of ambiguity parameters are ordered if V 2 (l/m(x)) 
has a definite sign. We prove these statements by using the Hellmann-Feynman theorem and 
offer examples of their application. 



1 Introduction 

Last few decades, quantum mechanical systems with position-dependent mass (PDM) 
have received considerable attention. The interest stems mainly from the relevance 
of the PDM background for describing the physics of compositionally graded crys- 
tals [HE] and semiconductor nano devices [3H5]. These applications have stimulated 
the study of the various theoretical aspects of the PDM Schrodinger equation; in par- 
ticular, its exact solvability [SHE], shape invariance [9], supersymmetry and intertwin- 
ing properties [T0H12] . point canonical transformation [HUH], iterative solution [15] . 
and relation to theories in curved spaces [16] have been examined. 

However, it is known that the PDM Schrodinger equation suffers from ambiguity 
in operator ordering, caused by the non-vanishing commutator of the momentum 
operator and the PDM. The PDM Hamiltonians with different ambiguity parameters 
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have been proposed [TTl - f2U] . but none of them can be preferred according to the 
existing reliability tests [2TH23] . Therefore the attempts are made to settle the issue 
by fitting the calculated binding energies to the experimental data [2^|l25]. 

For generelizing such findings and obtaining additional information, one needs 
some tools to compare the energy eigenvalues predicted by the different PDM Hamil- 
tonians. Within the constant-mass framework, a convenient tool is provided by the 
so-called comparison theorems [26-28J. For example, the elementary comparison the- 
orem [26l|28] states that if two real potentials are ordered, < V^ 2 \ then each 
corresponding pair of eigenvalues is ordered, < E^ 2 \ 

The purpose of this paper is to establish the comparison theorems that confront 
the energy eigenvalues of the constant-mass and PDM Schrodinger equations, as well 
as the energy eigenvalues of the PDM problems with different ambiguity parameters. 
Our presentation is based on the Hellmann-Feynman theorem [3U] and makes use of 
the ideas developed for the constant-mass case [281129] . 

The plan of the paper is as follows. In Section [2] we introduce the PDM Hamil- 
tonians and recall the Hellmann-Feynman theorem. In Section [3] the comparison 
theorems on the PDM background are formulated and proved. In Section H] we apply 
these theorems to two PDM problems of current interest. Finally, our conclusions are 
summarized in Section 

2 Preliminaries 

For the PDM Schrodinger equation the most general form of the Hamiltonian is given 
by HZ! 

H pdm = -i (m(x) a Vm(x) /3 Vm(x) 7 + m(x) 7 Vm(x) /3 Vm(x) Q ) + V(x) (1) 

where a, (3, 7 are the ambiguity parameters (a + /3 + 7 = —1) and the units with 
h = 1 are used. In this paper we will adopt the sets of the ambiguity parameter 
values suggested by BenDaniel and Duke [18] (a = 7 = 0, j3 — —1), Li and Kuhn [19] 
(a = P = -1/2, 7 = 0), and Gora and Williams [2DJ (a = -1, = 7 = 0). 

The methods we are going to apply are valid for arbitrary dimension N. We 
suppose that the Hamiltonian operators have domains T>(H) C L 2 (R N ), they are 
bounded below, essentially self-adjoint, and have at least one discrete eigenvalue at 
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the bottom of the spectrum. 

To derive our main results, we need the Hellmann-Feynman theorem [30]. This 
theorem states that if the Hamiltonian of a system is H(a), where a is a parameter, 
and the eigenvalue equation for a bound state is H(a)\a) = E(a)\a), where E(a) is 
the energy and \a) the normalized associated eigenstate, then 



dE(a) 
da 



dH{a) 



da 



Note that the proof relies on the self-adjointness of H(a) and does not change for 
PDM Hamiltonians. 

3 Comparison theorems 

First, let us formulate the theorem that confronts the energy eigenvalues of the 
constant-mass and BenDaniel-Duke PDM Hamiltonians with the same potentials. 

Theorem 1 

Suppose that the Hamiltonian 

# (0) = -tt^V 2 + V(x) (3) 
2m 

with a real potential V(x) and a constant-mass mo has discrete eigenvalues 
characterized by a set of quantum numbers {n}. Then the corresponding eigenvalues 
of the BenDaniel-Duke PDM Hamiltonian 

satisfy 

£{°}<£w D) if Vx 0<m(x)<m , (5) 
E^ } >E^ if Vx m(x)>m , (6) 

provided that these eigenvalues exist. 
Proof: Define the Hamiltonian 

H(a) = (1 -a)H {0) + a# (BD) , (7) 

which turns into and if ( BD ) when a = and a = 1, respectively. Assume that 
H(a) possesses well defined eigenvalues Ei n \(a), for < a < 1, and the normalized 
associated eigenfunctions in the coordinate representation are ■?/>{«} ( x > a )- 



Applying the Hellmann-Feynman theorem (J2J), we get 



/ (— U - — J |W> {n} (x;a)| 2 dx. (9) 
J V m ( x ) m oJ 



da 

where the integration is performed over the whole space and the asterisk denotes 

complex conjugation. 

Integrating by parts and taking into account that ip{ n \ (x; a) and Vip{ n } (x; a) must 

vanish at infinity, we obtain 

dE {n} (a) _ 1 ffl 1 \ 
da 2 

It is a positive (negative) number if < m(x) < m (m(x) > m ) for all x, so that 
Ei n \(a) is an increasing (decreasing) function of a. For definiteness, let < m(x) < 
m . Then it follows immediately 

^W(O) = E® < < } D) = E {n} (l) (10) 

that completes the proof. Note that an alternative proof can be given by applying 
the variational characterization j2I] of the discrete part of the Schrodinger spectrum. 

It is now tempting to compare the eigenvalues of the constant-mass Hamiltonian 
with those of PDM Hamiltonians other than the BenDaniel-Duke one. However, in 
that case we encounter an obstacle that becomes clear if we first find out how the 
eigenvalues of different PDM Hamiltonians are ordered. This is done in the following 
theorem. 



Theorem 2 

The discrete eigenvalues -E^^, and E^y of the BenDaniel-Duke, Li-Kuhn 

and Gora- Williams PDM Hamiltonians 

" ,BD, = -5 v ^ v + ,/M - < n > 

i/« LK » =--( -J = V^ = V + v—^v—!—) + v(x), (12) 
4 lym(x) A/wi(x) ym(x) ym(x) J 

H (GW) = _1( _1_ V 2 + V 1_L_) + r(x ) (13) 

satisfy 



4 \m(x) m(x) 



BS" < Bgf < BST 'f Vx VM-^)<0, (14) 



E^>E^>E^ if Vx V<(^)>0, 



m^x) 
1 



provided that these eigenvalues exist. 

Proof: Let us prove the inequalities for E^ and E^ . We define the parameter- 
dependent Hamiltonian H(a) by 

H(a) = (l-a) + a # (LK) (16) 



and make use of the Hellmann-Feynman theorem ([2]), to obtain 

dE {n} (a) 



=V + V- 



da 

1 _ I 



x; a) 



2 mix.) 



y/m(x) a/ m(x) ^/m(x) ^/m(x) 

Integration by parts yields 

-|W> {n} (x;a)| ax 



V>{ n }(x;a)o?x. (17) 



9a 2 / mfx 



4 J \V m ( x ) / v m ( x ) 

+ I / V f A-t ^Cx; a) J ■ J_ v ^fn}( x ; a)rfx 
4 y V Jm x / \/m(x) 



— i-ro (Vm(x)) • V|^{n}(x; a)| 2 dx 

' ^V 2 (— M |V> } (x;a)| 2 dx. (18) 

Let V 2 (l/m(x)) < for all x, then E{ n }(a) is an increasing function and we get 

E {n m = E^<E^f = E {n} (l). (19) 

that completes the proof. For the case of E^ and E^ , the proof is identical since 
the factor V 2 (l/m(x)) arises in this well. 

It is now evident from ([9]) and (|T8|) that if we try to compare e9^ with the 
constant-mass energy E^\, the sign of the integral will be determined by the signs 
of both (l/m(x) — l/m ) and V 2 (l/m(x)). Unfortunately, this leads to inconsistent 
conditions. For example, in order to get the inequality E]^\ > E^T 1 , we have to 
put l/m(x) < l/m and V 2 (l/m(x)) > 0, i.e., l/m(x) must be bounded from above 
and convex that is impossible. The same obstacle is encountered when dealing with 

E (GW) 
{n} 



5 



4 Applications 



In this section, we consider two specific PDM problems, which are discussed in litera- 
ture, and show how the comparison theorems explain the peculiarities of their energy 
spectra. 

4.1 Case 1 

The three-dimensional mass distribution of the form 

m W = n r° v» » ( 2 °) 

(1 + /err 



with r = |x| and nonnegative k, has been shown [16] to give rise to an exactly- 
solvable extension of the Coulomb problem, V(r) = —Ze 2 /r. This extension is useful 
as it enables one to trace the link between the PDM background and theories with 
deformations in the quantum canonical relations or with curvature of the underlying 
space. 

For this case, the discrete energy eigenvalues of the PDM Hamiltonian (JTJ, in units 
with h = m = e = 1, are written as [16] 



[Z-t(/(Z + l)-2/3)] 2 Zk 
E = 2^ + T" 



+ — [21(1 + 1) - n 2 - 4/3 + (1 + 4a)(l + 4 7 )1 (21) 

o 

where I = 0,1, ... and n = l+l, 1+2, ... are the orbital and principal quantum numbers, 
respectively. In contrast to the constant-mass Coulomb problem, the system has only 
a finite number of discrete levels, so that the allowed values of I and n are restricted 
by 

^[l(l + l) + n 2 -2(3} <Z. (22) 

Such a restriction implies that in presence of the PDM the energy eigenvalues may 
be closer to continuum and thus larger than the ordinary Coulomb eigenenergies 
= -Z 2 /{2n 2 ) calculated with the mass m . 
It is Theorems 1 and 2 that permit us to determine how the energy eigenvalues are 
ordered. Since in f[2"U|) we have m(r) < m , the eigenvalues of the BenDaniel-Duke 
PDM Hamiltonian must obey _E( BD ) > E^°\ by Theorem 1. Since V 2 (l/m(r)) = 
(/t/mo)(6/€ + 4/r) > 0, it follows from Theorem 2 that the eigenvalues of the PDM 
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Hamiltonians with different ambiguity parameters are ordered as E^ BB ^ > E^ hK ' > 

£(GW)_ 

In order to illustrate these inequalities, we present figured] where we plot the energy 
for the ground state (n = 1, / = 0) and the first radially excited state (n = 2, 1 = 0), 
as a function of the deforming parameter k. In figure [I] the solid lines correspond 
to the constant-mass case whereas the broken curves represent the PDM cases with 
different ambiguity parameters. The circles indicate the points at which the bound 
states disappear according to (122]) . From figure [I] we see that, for all allowed k, it 
holds E^ BD ^ > E(°\ as it was proved, and also E^ LK ^ > E^°\ but we observe both 

E (GW) > E (0) and £(GW) < E (0) 

regions. Furthermore, we can see that the second 
proved inequality, £ (BD) > £ (LK) > E^ GW \ is indeed fulfilled. 
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Figure 1: Energy for the states (a) (n = 1, I = 0) and (b) (n = 2, I = 0), in the Coulomb 
potential and the mass distribution ([20]) . calculated with Ti = niQ = Ze 2 = 1. The solid line is the 
constant-mass result, the dotted, dashed and dash-dotted curves are the PDM results obtained with 
Hamiltonians (11), (12) and (13), respectively. 



(b) 



4.2 Case 2 

Now let us consider the one-dimensional mass distribution 

m(x) = m (1 + kx 2 ), (23) 

which is found to be useful for studying quantum wells [3]. Applying Theorem 1 to 
this PDM profile, we get the inequality _E( BD ) < E^ that justifies the shift of electron 
and hole binding energies to lower values which was observed in [3] when the spatial 
dependence of mass was included. On the other hand, Theorem 2 does not apply 
since the quantity V 2 (l/m(x)) = k,(6kx 2 — 2)/m (l + kx 2 ) 3 has an indefinite sign. 
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It is worth examining how this sign indefiniteness affects the energy spectrum. To 
that end, we choose the harmonic-oscillator potential, V(x) = ^mou 2 x 2 , for which 
the accurate numerical solution of the PDM Schrodinger equation with the mass 
distribution (123]) is available [15] . In figure |2] we plot the corresponding energy of the 
ground and the fifth excited states, as a function of k, for the three PDM Hamiltonians 
with different ambiguity parameters. The energies have been calculated with U = 
m — oj — 1, by using the shooting method, and are in agreement with those computed 
in [15] where the results obtained with the same mo and u, and k — 0.1 are reported. 
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Figure 2: Energy for the states (a) n = and (b) n = 5, in the harmonic-oscillator potential and 
the mass distribution ([2"3")l . calculated with Ti = mo = oj = 1. The dotted, dashed and dash-dotted 
curves are the PDM results obtained with Hamiltonians (11), (12) and (13), respectively. 



From figure [2] it is evident that for the excited state the discrepancy among the 
energies evaluated using the different PDM Hamiltonians is less profound. However, 
we call attention to a serious difference between the ground and excited states. As 
seen in figure [2J the ground-state energies are ordered as E^ BD ^ < £"( LK ) < E^ GW ^ 
whereas the energies of the fifth excited state (and of the states with n > 5) are in 
inverse order. This inversion can be easily understood in conjunction with Theorem 
2. It is known that the wave functions of highly excited states are spread to larger 
distances. Consequently, with increasing n, the mean value of x 2 grows and eventually 
reaches the point where the sign of V 2 (l/m(x)) in Theorem 2 reverses, thus inverting 
the order of energies. 



8 



5 Summary 

In this paper, we have established the comparison theorems for the PDM Schrodinger 
equation. Our first theorem states that the corresponding eigenvalues of a constant- 
mass Hamiltonian and of a BenDaniel-Duke PDM Hamiltonian with the same poten- 
tial are ordered if the constant and position-dependent masses are ordered everywhere. 
The second theorem concerns PDM Hamiltonians with the different sets of ambiguity 
parameters: the BenDaniel-Duke, Li-Kuhn, and Gora- Williams Hamiltonians. It is 
proved that their corresponding eigenvalues are ordered if the Laplacian of the inverse 
mass distribution l/m(x) has a definite sign. 

We have applied these theorems to the PDM Coulomb and harmonic-oscillator 
problems and have been led to the following conclusions. First, the eigenvalues of 
PDM Hamiltonians other than the BenDaniel-Duke one do not have to be in the 
strict order with respect to the eigenvalues of the constant-mass Hamiltonian. For 
instance, from both figures [T] and [2] it is seen that the order of the Gora- Williams and 
constant-mass ground-state energies do vary, depending on the value of the deform- 
ing parameter k. Second, if the quantity V 2 (l/m(x)) has no definite sign and thus 
Theorem 2 does not apply, the order of the energies calculated using different PDM 
Hamiltonians may alternate, as seen by comparing parts (a) and (b) of figure [2j We 
therefore think that for establishing further comparison rules within the PDM frame- 
work one should restrict the potential profile to, e.g., a spherically-symmetric case, 
the way the generalized comparison theorems for the ordinary Schrodinger equation 
have been obtained [2"T] . 

The comparison rules we have found out can be employed for analyzing the energy 
spectra in semiconductor nano devices; an example of application to the quantum 
well system was sketched in the previous section. In this connection, it is worthwhile 
to extend the present approach to periodic heterostructures, which allow the direct 
fit of PDM binding energies to experiment [25] • Then we will have to abandon the 
requirement of vanishing of the wave function at infinity which the proof of our 
theorems relies on. What comparison rules might be formulated in that case is an 
interesting open question. 
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